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Abstract. We present a study of the dynamics and magnetic activity of M dwarfs using the 
largest spectroscopic sample of low-mass stars ever assembled. The age at which strong surface 
magnetic activity (as traced by Ha) ceases in M dwarfs has been inferred to have a strong 
dependence on mass (spectral type, surface temperature) and explains previous results showing 
a large increase in the fraction of active stars at later spectral types. Using spectral observations 
of more than 40000 M dwarfs from the Sloan Digital Sky Survey, we show that the fraction of 
active stars decreases as a function of vertical distance from the Galactic plane (a statistical 
proxy for age), and that the magnitude of this decrease changes significantly for different M 
spectral types. Adopting a simple dynamical model for thin disk vertical heating, we assign an 
age for the activity decline at each spectral type, and thus determine the activity lifetimes for 
M dwarfs. In addition, we derive a statistical age-activity relation for each spectral type using 
the dynamical model, the vertical distance from the Plane and the Ha emission line luminosity 
of each star (the latter of which also decreases with vertical height above the Galactic plane). 

Keywords, stars: activity, stars: kinematics, stars: late-type,stars: low-mass, brown dwarfs 



1. Introduction 

M dwarfs are the smallest and least luminous stars on the main sequence, yet they 
are the most numerous of all stellar constituents and have lifetimes longer than the 
current age of the Universe. Determining the age of an individual main sequence, field M 
dwarf is quite challenging. Stellar clusters would seem to provide the ideal environments 
for calibrating the ages of M dwarfs (because the age of the cluster is determined). 
However, the large distances to clusters older than a few Gyr (and resulting faintness 
of cluster M dwarfs), preclude any detailed observations of the intermediate or oldest 
cluster resident M dwarfs. Age estimates for field M dwarfs must therefore rely on a 
number of techniques that piggyback on other age dating methods. Although several 
methods exist for estimating the ages of M dwarfs (e.g. contribution by S. Catalan), we 
focus this contribution on what the magnetic activity of M dwarfs can tell us about their 
age. 

M dwarfs are host to intense magnetics dynamos that give rise to chromospheric and 
coronal heating, producing emission from the X-ray to the radio. This magnetic activity 
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has long been thought to be tied to the age of the host star. Almost 40 years ago, Wilson 
& Woolley (1970) found a link between magnetic activity (as traced by the Ca II H & K 
emission lines) and the orbital elements (namely the inclination and eccentricity) of more 
than 300 late-type dwarfs; stars in near-circular orbits with small inclinations had the 
strongest activity. They concluded that as stars age, their orbits get more inclined and 
more eccentric (due to dynamical interactions), and thus, magnetic activity must also 
decline with age. Subsequent studies over the following decades have found similar con- 
nections between age and activity in low-mass stars (Wielen 1977; Giampapa & Liebert 
1986; Sodcrblom et al. 1991; Hawley ct al. 1996, Hawley et al. 1999, Hawley ct al. 2000). 

Although activity has been observed for decades, the exact mechanism that gives rise 
to the chromospheric heating is still not well-understood. In the Sun, magnetic field gen- 
eration and subsequent heating is closely tied to the Sun's rotation. From hclioscismology, 
we know that there is a rotational boundary between the inner solid-body rotating radia- 
tive zone and the outer differentially rotating convective zone (Parker 1993; Ossendrijver 
2003; Thompson et al. 2003). This boundary, dubbed the tachocline, creates a rotational 
shear that likely allows magnetic fields to be generated, preserved and eventually rise to 
the surface where they emerge as magnetic loops. These loops bring heat to the Sun's 
chromosphere and corona, driving both large stellar flares, as well as lower- level quies- 
cent (or persistent) magnetic activity. The faster a star (with a tachocline) rotates, the 
stronger its magnetic heating and surface activity. Therefore, the angular momentum 
evolution of solar-type stars should play an important role in determining the magnitude 
of the observed activity. 

Angular momentum loss from magnetized winds has been shown to slow rotation in 
solar-type stars; as a result magnetic activity decreases. Skumanich (1972) foimd that 
both activity and rotation decrease over time as a power law (t^'^'^). Subsequent studies 
confirmed the Skumanich results and empirically demonstrated a strong link between age, 
rotation and activity in solar-type stars (Barry 1988; Sodcrblom ct al. 1991; Pizzolato et 
al. 2003; Mamajek & Hillenbrand 2008; sec contribution by E. Mamajek). 

There is strong evidence that the rotation-activity (and presumably age) relation ex- 
tends from stars more massive than the Sim to smaller dwarfs (Pizzolato et al. 2003; 
Mohanty & Basri 2003; Kiraga & Stepien 2007). Therefore, rotation derived ages (using 
gyrochronology) may provide a useful independent estimate of age for M dwarfs when 
large enough calibration samples can be acquired (see contributions by S. Meibom, J. Ir- 
win, and S. Barnes). However, at a spectral type of ^M3 (0.35 Mq; Reid & Hawley 2005; 
Chabrier & Baraffe 1997), stars become fully convective and the tachocline presmnably 
disappears. Despite this changes, magnetic activity persists in late-type M dwarfs; the 
fraction of active M dwarfs peaks around a spectral type of M7 before decreasing into 
the brown dwarf regime (Hawley et al. 1996; Gizis et al. 2000; West et al. 2004). 

It is unclear if the rotation-age-activity relation extends to the fully convective M 
dwarfs. A few empirical studies have uncovered evidence that activity and rotation might 
be linked in late-type M dwarfs (Delfosse et al. 1998; Mohanty & Basri 2003; Reiners & 
Basri 2007). In addition, recent simulations of fully convective stars find that rotation may 
play a role in magnetic dynamo generation (Dobler et al. 2006; Browning 2008). However, 
a recent study found detectable rotation in a few inactive M7 dwarfs, indicating that the 
situation may in fact be more complicated (West & Basri 2009). 

Irrespective of rotation, many studies have found evidence that the age-activity relation 
extends into the M dwarf regime. Eggen (199) observed a Skumanic;li-typc decrease in 
activity as a function of age. Larger samples of M dwarfs have added further evidence that 
magnetic activity decreases with age (Fleming et al. 1995; Gizis et al. 2002). There is also 
evidence that M dwarfs may have finite active lifetimes. Stauffer et al. (1994) suggested 



Activity Ages of M Dwarfs 



121 



that activity may not be present in the most massive M dwarfs in the Pleiades. Hawley 
et al. (2000) confirmed the Stauffer et aL (1994) claim by observing a a sample of clusters 
that spanned several Gyr in age. They were able to calibrate the lifetimes for early-type 
M dwarfs by observing the color at which activity (traced by Ha emission) was no longer 
present. Because of the small sample size, the derived Hawley et al. (2000) activity ages 
are lower limits at a given color or spectral type. As mentioned above, the Hawley et al. 
(2000) study was limited to younger, nearby clusters due to the intrinsic faintness of M 
dwarfs; thus, it could only probe ages of a few Gyr and spectral types as late as ^^MS. 
Larger samples of M dwarfs are required to statistically derive age-activity relations that 
extend both in the ages they probe and the range of spectral types they cover. 

In this contribution, we review recent studies that have used large spectroscopic sam- 
ples to investigate the relationship between age and activity in M dwarfs. In addition, 
we derive a statistical Ha age- activity relation for M2-M7 dwarfs. 

2. SDSS DR5 Low-Mass Star Spectroscopic Sample 

Large surveys such as the Sloan Digital Sky Survey (SDSS; York et al. 2000) have 
created optical and infrared catalogs of tens of millions of M dwarfs (Bochanski et al. 
2009 in prep). In addition to the photometric data, the SDSS has obtained spectra for 
tens of thousands of M dwarfs. Recently, these surveys have been utilized to examine the 
statistical properties of M dwarfs, the dynamics of the Milky Way and detailed studies 
of magnetic activity (Hawley et al. 2002; West et al. 2004; West et al. 2006; Covey et al. 
2007; Bochanski et al. 2007a, 2007b; West et al. 2008; hereafter W08). 

Recently, the SDSS Data Release 5 (DR5) low-mass star spectroscopic sample was 
released (wosjil, containing over 44000 M and early L-type dwarfs. Radial velocities 
(accurate to within '-~^5 km s^^) and photometric distances (from the r — z colors; see 
Bochanski et al. 2009 in prep) were measured for all of the stars. Proper motions (matches 
to USNO-B; Munn et al. 2004) were obtained for over 27000 of the stars, allowing for full 
U, V, W 3-D space motions for over 27000 low-mass dwarfs, the largest such sample ever 
assembled. In addition, all of the spectra were run through the Hammer spectral typing 
facility (Covey et al. 2007) to assign spectral types and measure the strength of the Ha 
emission line. One of the hallmarks of the W08 sample is that it probes the entirety of the 
Galactic thin disk and extends well into the thick disk (~50-2500 pc above the Plane). 
This is due to the fact that the majority of SDSS fields are in the North Galactic Cap 
and that SDSS was designed to study the distant Universe; intrinsically faint Galactic 
M dwarfs are bright in a deep extragalactic survey. 

3. Results 

3.1. Ha Activity 

Activity in M dwarfs varies as a function of spectral type. Figure [1] shows the fraction 
of active M dwarfs as a function of spectral type from West et al. (2004). The active 
fraction is very small for early- type M dwarfs, peaks at M7-M8 and declines into the L 
dwarf regime. While some of the morphology of this relation may be due to the ability 
(or inability) to host a strong dynamo (the activity fraction decrease in the late-type M 
and L dwarfs is likely caused by the atmospheres becoming neutral), a large part of the 
shape may be due to age effects. If activity has a finite lifetime that changes as a function 

t Measured quantities can be obtai ned electronically using the CDS Vizier database 
http://vizier.u.strasbg.fr/vis-bin/VizieR| 
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of spectral type, and we observe a range of stellar ages, then stars with longer lifetimes 
will appear to have higher activity fractions. 
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Figure 1. Fraction of active stars as a function of spectral type (reproduced from West et al. 
2004). Numbers above each point represent the number of stars used to compute the fraction. 

To test this hypothesis. West et al. (2004, 2006) showed that the fraction of active M7 
stars decreases as a function of vertical distance from the Galactic Plane (see Figure [2]). 
Stars are born in dynamically cold molecular gas near the midplane of the Galaxy. Over 
time they undergo dynamical interactions with molecular clouds, which add energy to the 
stellar orbits in all dimensions (see contribution by B. Nordstrom). It is this process that 
gives thickness to the Milky Way disk and allows us to use vertical distance as a proxy 
for age; stars further from the Plane are statistically older (they have to be dynamically 
heated to reach those heights), while stars near the Plane are statistically younger. Figure 
[2] shows that the younger stars, near the plane of the Milky Way are almost all active, 
whereas the fraction of active stars at larger distances above (and below) the Galaxy is 
significantly reduced; older stars have ceased being active. The decrease of active fractions 
as a function of vertical height above the Plane is seen for all M dwarf spectral types 
(W08). 
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Figure 2. The fraction of active M7 stars as a function of vertical distance from the Galactic 
plane (reproduced from West et al. 2006). There is a significant decrease in the active fraction 
as a function of Galactic height, which can be used as a proxy for age. Younger stars near the 
Plane are almost all active, whereas older stars, further from the Plane have ceased being active. 
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W08 also showed that the amount of activity, quantified by LnaZ-^boi (the ratio of 
the luminosity in the Ha emission line to the bolometric flux of the star) decreases over 
time. Figure [3] shows the median log(LHa/^boi) as a function of vertical distance from the 
Plane. The narrow error bars represent the spread of the values and the wide error bars 
indicate the uncertainty of the median relation in each bin. The decrease as a function of 
height is a statistically significant over most of the spectral types. Some of the spectral 
types were omitted from Figure [3] because they lacked a sufficient number of active stars 
for a robust study. 
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Figure 3. Median log(LHa/-i'boi) as a function of vertical distance from the Plane (reproduced 
from W08). The narrow error bars represent the spread of the values and the wide error bars 
indicate the uncertainty of the median relation in each bin. The decrease as a function of height 
is a statistically significant over most of the spectral types. 



3.2. 1-D Dynamical Model 

West et al. (2006, 2008) devolved a ID dynamical model that traces the vertical dynamics 
of stars as a function of time. The model assumes a constant star formation rate and 
adds a new population of 50 stars every 200 Myr for a total simulation time of 10 Gyr. 
Each new population of stars begins with a randomly drawn velocity dispersion 8 km 
s~^ (Binney et al. 2000). Orbits are integrated using a "leap frog" integration technique 
(Press et al. 1992) and the vertical Galactic potential from Siebert et al. (2003). 

Dynamical heating was simulated by altering the velocities of stars such that their 
velocity dispersions as a function of age match a aw oc relation (Wielen 1977; Fuchs 
et al. 2001; Hanninen & Flynn 2002; Nordstrom et al. 2004). Energy was added to stars 
that were within a certain distance from the Plane. This "region of infiuence" is a way 
to parameterize the cross section of interaction with molecular clouds during a Plane 
crossing (see West et al. 2006 for more information). The majority of the molecular gas 
is constrained to a small range of Galactic heights and the dynamical interaction depends 
on the proximity to the cloud, the density of the gas and the velocity of the star, all of 
which are absorbed in the "region of influence" . The size of the "region of influence" was 
varied from ±0.5 to ±5 pc in intervals of 0.5 pc. Seventy simulations were run for each 
spectral type. Each simulation recorded the velocity, position, and age of each star. 

3.3. Age-Activity 

W08 introduced a binary activity state to the dynamical models; stars started their lives 
as active and after a finite "lifetime" , the activity turned off. The active lifetimes were 
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varied from 0.0 to 9.0 Gyr in 0.5 Gyr intervals for each spectral type. The resulting model 
active fractions (as a function of vertical distance from the Plane) were compared to the 
empirical relations using a chi-squared minimization technique and a best-fit model for 
each spectral type was determined. The resulting Ha activity lifetimes are shown in 
Figure S) The results from the Hawley et al. (2000) cluster study are overplotted for 
comparison (dotted). W08 found that there is a significant increase in activity lifetimes 
between spectral types M3 and M5, possibly indicating a physical change in the produc- 
tion of magnetic fields as the stellar interiors become fully convective. 
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Figure 4. The Ha activity lifetimes of M dwarfs as determined by comparing the SDSS DR5 
spectroscopic data to ID dynamical simulations (stars; reproduced from W08). The Hawley 
et al. (2000) activity lifetime relation is overplotted for comparison. As predicted, the Hawley 
relation provides a lower limit to the ages. W08 found that there is a significant increase in 
activity lifetimes between spectral types M3 and M5, possibly indicating a physical change in 
the production of magnetic fields as full convection sets in. 

Figure [3] suggests that activity is not simply a binary process, but rather that the 
amount of magnetic activity may decrease over time. To quantify this decrease, we com- 
pared the median log(LHa/iboi) from Figure [3] to the median ages in the same vertical 
distance bins from the W08 ID dynamical models. The result, shown in Figure [5l gives 
a statistical, quantitative relationship between the Ha activity of an M dwarf and it's 
age. Horizontal error bars represent the spread in age at a given bin. The age-activity 
relations are consistent with a smooth Skumanich-like activity decrease until the star 
reaches its activity lifetime, after which it rapidly falls to its eternal inactive state. 

We fit the activity-age relations in Figure [5] with a function of the form: 

log(LHa/^boi) = J^7—^ - ^' (3-1) 

where a, b and n are coefficients, I is the active lifetime (Gyr), and t is the age measured 
in Gyr. The exponent n was forced to be the same for all spectral types. The lifetime 
/ coefficient was allowed to vary within the uncertainties of the derived active lifetimes 
from W08 (Figure [3]). The resulting best-fit coefficients can be found in Table [TJ We 
caution that these relations are purely statistical and may not be appropriate for use 
with individual stars. Additional discretion should be used with the functional fits, which 
may not be justifiably extrapolated beyond where data exist. 
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Figure 5. Median log(LHa/-I/boi) as a function of age (as derived from the median age of the same 
bins in the W08 dynamical models). The error bars represent the spread in the distributions. 
The best-fit function is overplotted for each spectral type (dashed). 



Table 1. Coefficients for best-fit activity-age relation 



Spectral Type 


a 


b 


n 


' (Gyr) 


M2 


0.106 


3.38 


2.0 


1.8 


M3 


0.213 


3.54 


2.0 


2.7 


M4 


1.41 


3.39 


2.0 


5.0 


M5 


2.85 


3.57 


2.0 


7.0 


M6 


1.78 


3.90 


2.0 


7.0 


M7 


11.8 


4.01 


2.0 


8.0 



4. Summary 

Large astronomical surveys have produced M dwarf samples of unprecedented size. 

We have shown that using the statistical foothold of the largest low-mass spectroscopic 
sample ever assembled, a strong tie between magnetic activity (as traced by Ha) and 
age has been established and that M dwarfs have finite active lifetimes. These lifetimes 
are a strong function of spectral type, dramatically increasing as the M dwarf interiors 
become fully convective. We derived age-activity relations for M2-M7 dwarfs based on 
the average activity as a function of height above (or below) the Galactic plane. Future 
studies will extend these relations to other spectral types. In addition, these relations can 
be confirmed and calibrated using both white dwarf cooling ages in wide binary systems 
that have both a white dwarf and an M dwarf (sec contributions by S. Catalan, M. 
Salaris, and J. Kalirai; Silvcstri ct al. 2006), and deep cluster observations of M dwarfs 
(when such observations become possible). 
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